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Abstract. An introduction to nucleosynthesis, the creation of the elements 
in the big bang, in interstellar matter and in stars is given. The two-step pro- 
cess *He(2n,7)®He and the reverse photodisintegration ^He(7,2n)*He involv- 
ing the halo nucleus ®He could be of importance in the a-process in type-II 
supernovae. The reaction rates for the above processes are calculated using 
three-body methods and show an enhancement of more than three orders of 
magnitude compared to the previous adopted value. Direct-capture calculations 
give similar values for the above reaction rates. Therefore, this method was also 
used to calculate the reaction rates of the two-step processes '^He(2n,7)*He and 
®Li(2n,7)^'^Li and the reverse photodisintegration of *He and ^^Li that could 
be also of importance in the a-process. 



1 Introduction 

One of the main driving forces for the evolution of our universe are nuclear 
processes. These nuclear processes are responsible for the energy production 
in stars as well as for the production of the elements in our universe. The 
process of the creation of the elements is called nucleosynthesis. We have three 
different scenarios, where nucleosynthesis occurs: in the big bang (primordial 
nucleosynthesis) , in interstellar space (interstellar nucleosynthesis) and in stars 
(stellar nucleosynthesis). The result of nucleosynthesis is reflected today by 
the observed elemental abundances in different objects (sun and other stars, 
planets, meteorites, interstellar matter, . . . ). 

The light elements from hydrogen to lithium with mass numbers A — 1 to 
A = 7 were mainly synthesized in primordial nucleosynthesis. In interstellar 
nucleosynthesis primarily elements from lithium to boron with mass numbers 
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A = 6 to A = 11 were produced. Finally, all the heavier elements from car- 
bon to uranium (in the astrophysical jargon these elements are called metals) 
with mass numbers A = 12 to A = 238 were created in stellar nucleosynthe- 
sis. Even though the heavier elements amount to about only 1% of the total 
observed abundance they are obviously indispensable for our environmemt. 
Without stellar nucleosynthesis our universe would be a boring place with the 
two gases hydrogen and helium (with a negligible amount of lithium, beryllium 
and boron) and without all the elements necessary for the existence of life. 

One of the essential inputs for the investigation of different astrophysical 
scenarios for nucleosynthesis are cross sections, astrophysical S~factors and re- 
action rates of nuclear reactions as well as half lives and decay rates of nuclear 
decays. These quantities are determined using different experimental facilities 
(accelerators, reactors, . . . ). However, in many cases it is not possible to obtain 
the needed cross sections directly from experiments. This is due to the fact that 
astrophysically relevant reactions take place mainly at thermonuclear energies. 
These energies are mostly well below the Coulomb and/or centrifugal barrier, 
and in this region the cross sections are often too small for experimental deter- 
mination. Furthermore, in nucleosynthesis often radioactive nuclei are involved, 
for which experimental information is only scarcely available. For nuclei far-off 
stability often experimental information is not attainable at all. 

In the next section a short introduction to nucleosynthesis is given em- 
phasizing those astrophysical scenarios, where light nuclei are involved and 
methods developed in the few-body field could be applied. In Sect. 3 we show 
as an example for an application of few-body calculations for astrophysically 
relevant nuclear reactions the two-step process '*He(2n,7)^He. This reaction 
involves also the halo nucleus ^He. In the last section we calculate the reaction 
rate for this reaction as well as other two-step reactions involving halo nuclei, 
like ^He(2n,7)^He and ^Li(2n,7)-'^^Li, using the simpler direct-capture (DC) 
model. 



2 Nucleosynthesis 

In this section we only can give a short introduction to nucleosynthesis. There 
are many textbooks where this subject can be found in more detail (e.g. |^). 

Primordial nucleosynthesis took place as the cooling early universe reached a 
temperature of about 10^ K about three minutes after the big bang. Before that 
time the temperature was high enough that formed deuterons were immediately 
again destroyed through photodisintegration. From primordial nucleosynthesis 
calaculations we can deduce abundance ratios of the light stable nuclides ^H, 
^H, ^He, ^He, ^Li and ^Li. Even some primordial abundance ratios may be 
altered somewhat by interstellar and stellar nucleosynthesis, the agreement 
between the abundances obtained from primordial nucleosynthesis with the 
observed abundances is remarkable. 

In interstellar nucleosynthesis the nuclides ^Li, ^Li, ^Be, ^°B and ^^B are 
produced. In this scenario high-energy (larger than about 1 GeV) cosmic rays 
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hit interstellar matter and create through spallation processes the above nuclei. 
Table 2.1. Nuclear burning phases and corresponding temperatures 



Nuclear burning phase 


Temperature (Tg) 


Hydrogen burning 


0.01-0.04 


Helium burning 


0.1-0.2 


Carbon burning 


0.6-0.8 


Neon burning 


1.2-1.4 


Oxygen burning 


1.5-2.2 


Silicon burning 


3-4 



Stellar nucleosynthesis begins when through gravitational contraction in a 
star the density and temperature in the core gets high enough so that nuclear 
processes are ignited. Stellar nucleosynthesis takes place in succesive burning 
phases. When the fuel of a preceding burning phase has been consumed, the ash 
of this burning phase is the fuel for the next burning phase. In the outer layers 
the previous burning phases still continue. The different burning phases are 
hydrogen burning (pp-chain, CNO-cycle, Ne-Na cycle, Mg-Al cycle), helium, 
carbon, neon, oxygen and silicon burning. In hydrogen burning helium is pro- 
duced, in helium burning carbon and oxygen are synthesized, in the remainig 
advanced burning phases elements up to elements in the region of iron are pro- 
duced. In successive burning phases the temperature must get larger, because 
of the increasing Coulomb barriers of the fusing heavier nuclei (Table ^.l[ ). 

The stars in which hydrogen burning takes place are named main-sequence 
stars, whereas the other burning phases occur in red giants. Not in all stars all 
the burning phases are ignited. For stars with less than approximately 8 solar 
masses, after helium burning, the outer part of the star that is blown off by 
strong stellar winds is called a planetary nebula, whereas the core remains as a 
white dwarf. For stars with more than about 8 solar masses all the above cited 
burning phases take place. Finally, after silicon burning the core of the star 
consists of elements with mass numbers in the region of iron. Then it is not 
possible to gain energy anymore through nuclear fusion and the star collapses. 
Through the rebounce effect after the implosion the stellar mantle is ejected in 
a gigantic explosion, whereas the core remains as a neutron star or black hole. 
This dramatic event is called a supernova of type II (a supernova of type I 
takes place in a binary system consisting of two stars, where triggered through 
mass accretion from the accompaning star the star is disrupted completely). 

There are astrophysical scenarios in which free neutrons play a dominat- 
ing role. Most elements heavier than iron have been created through neutron 
capture in two processes: in the s~process occuring in helium burning of red 
giants and in the r-proccess occuring in a supernova of type II. The neutron 
flux in the s-process of helium burning is comparetively low (iV,j ~ 10^ cm^"^), 
whereas in the r-process it is much higher {Nn — 10^°cm~'^). 
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Figure 1. Primordial nucleosynthesis in the Standard Big Bang (SBB) 



In Fig. |l| the nuclear reaction network for primordial nucleosynthesis in the 
Standard Big Bang (SBB) is shown. Fig. || shows the pp-chain dominating 
in low-mass main-sequence stars, like our sun, whereas the CNO-cycle takes 
over for more massive stars. There are two other scenarios additionally to the 
before discussed s- and r-process where neutron capture plays also a domi- 
nating role. The first one is an alternative to the standard Big Bang and is 
called the Inhomogcnous Big Bang (IBB) with neutron densities up to about 
TVn ~ 10^° cm^"^. The second one is the so-called a-process occuring in type-II 
supernovae with neutron densities of about — 10^°~^° cm~'^. Fig. ^ shows 
the net-work up to a mass number of ^ = 12 which should be taken into ac- 
count additionally to Fig. ^ for primordial nucleosynthesis in the IBB and the 
a-process. These are neutron-rich scenarios, where reactions involving light 



5 



p(p,e+v)d 



d(p,y)^He 




Chain 1 Chain 2 Chain 3 



Figure 2. The pp-chain in liydrogen burning of main-sequence stars 



nuclei on the neutron-rich side of the region of stabihty have to be included. 
The main nuclear reactions up to a mass number of A — 12 taking place in the 
SBB, IBB, pp-chain, helium burning and the a-process are listed in Table 

We want to discuss the a-process a little more, because the reaction rates 
calculated in the next section occur in this scenario. The a-process takes place 
in the neutrino-heated hot bubble between the nascent neutron star and the 
overlying stellar mantle of a type-II supernova |^, @, |^, ^ • During the grav- 
itational collapse in the supernova the temperature in the core of the star gets 
so high (up to 10^° K) that all the nuclei are disassembled through photodis- 
sociation into protons and neutrons. Afterwards, the material being originally 
in nuclear statistical equilibrium (NSE) at high temperature, is expanded and 
cooled so rapidly that not all the a-particles have time to reassemble. In this 
environment most nucleons are either in the form of free neutrons or bound in 
a-particles. The nucleosynthesis in the neutrino bubble takes place in subse- 
quent steps: first, a-particles are formed from the free nucleons. Then, in the 
following a-process, nuclei up to about A « 100 are produced ||]. Finally, the 
neutrino bubble is also an ideal site for the r-process synthesizing the elements 
of about A > 100 |, 0, |. 
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Table 2.2. Nuclear reactions of astrophysical interest for light nuclei up to ^ = 12. 
The left column lists the major reactions for nucleosynthesis in the Standard Big 
Bang (SBB), additional reactions in the pp-chain and helium burning. The right 
column lists reactions which can be additionaly activated in the a-process or the 
Inhomogenous Big Bang (IBB). 



Target 


Residual 


Type of 


Target 


Residual 


Type of 


Nucleus 


Nucleus 


Reaction 


Nucleus 


Nucleus 


Reaction 




SBB 






a-process, IBB 


In 




13- 


4He 


6He 


(2n,7) 


iH 




(n,7) 


^He 


8He 


(2n,7) 


2H 


3H 


(n,7), (d,p) 


^He 


6Li 


/3- 


2H 


^He 


(P,7), (d,n) 




lOg 


(a,7) 


2H 




(a. 7) 




8Li 


(n,7) 


3H 


3He 




^Li 


"B 


(a, 7) 


3H 


4He 


(d,n), (p,7) 


^Be 


SBe 


(n,7) 


3H 




(a, 7) 


8He 


8Li 


/3- 


3He 




(n,p) 


8Li 


9Li 


(n,7) 


3He 


4He 


(d,p), (n,7) 


8Li 


8Be 


/3- 


^Hc 


^Be 


(a, 7) 


8Li 


"B 


(a,p) 


6Li 




(n, a) 


«Bc 


9Bc 


(n,7) 


6Li 


3He 


(P,a) 


»Li 


9Be 


/3- 


6Li 


^Li 


(n,7) 


9Li 


"Li 


(2n,7) 


6Li 


^Be 


(P:7) 


9Bc 


6Li 


(P,a) 


^Li 


4He 


(P,a) 


»Bc 


lOBe 


(n,7) 


^Be 


^Li 


(e-,z/) 


9Be 


10 B 


(P,7), (d,n) 


^Be 


4He 


(n,Q!) 


9Be 


12b 


(a,p) 


^Be 


^Li 


(n,p) 


9Be 


12c 


(a, n) 




pp-chain 


lOBe 


"Be 


(n,7) 


iH 




(p,e+i/) 


lOBe 


lOg 


/?- 




3He 


(P>7) 


lOg 


^Li 


(n,a) 


3Hc 


^He 


(3He,2p) 


lOg 


"B 


(n,7) 


3He 


^Be 


(a. 7) 


"Li 


"Be 


/3- 


^Be 


SB 


(^,7) 


"Be 


i^Be 


(n,7) 




«Be 


/3+ 


"Be 


"B 


/3- 


8Be 


^Hc* 


decay 


"B 


12b 


(n,7) 




helium burning 


"B 


12c 


(P,7), (d,n) 


^He 


12c 


(2a, 7) 


i2Be 


12b 


/3- 


8Be 


4He 


decay 


123 


12c 


r, (P,n) 
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3 Few Body Calculations 

Few-body methods have been used only partially for the calulation of astropy- 
sical S-factors or reaction rates of astrophysically relevant nuclear reactions. 
An example is an exact three-nucleon calculation for radiative capture by 
deuterons starting from realistic NN-potentials. In this calculation the astro- 
physically relevant reaction rates for ^H(p,7)'^He at thermonuclear energies and 
^H(n,7)^H at thermal energies have been determined Another example is 
the calculation of the non-resonant part of the cross section for ^He(d,7)^Li in 
the thermonuclear energy region. It was calculated in an a+p+n model by solv- 
ing the three-body Faddeev equations ||l3 and using a three-body variational 



approach in the asymptotic region 1 1 1 



In the following we present the two-step processes '*He(2n,7)^IIe and the 
photodisintegration of ^He as an example for an application of few-body cal- 
culations to astrophysically relevant nuclear reactions. A more extensive pre- 
sentation of this calculation is given in ref. 12 . 

The reaction rate per particle triplet for^IIe(2n,7)^He, in analogy to the 
triple-alpha process, is given by 

(2n''He) = 2/ dEi — , 

\ / Jo 'r(5Hc,^i) dEi 

d/n^Hc) (Ei,E2) 
dE2 ' , 1 

dE2 ^ ' 

where Ei and E2 are the relative energies in the (*He-|-n)- and (^He-|-n)- 
system, respectively. To determine the reaction rate in Eq. (|^), an inte- 
gration over both energies Ei and E2 has to be performed. The quantity 
r (^He, El) is the energy-dependent width of ^He, whereas d (n^He) {Ei) /dEi 
and d (n^He) {Ei, E2) /dE2 arc the integrands of the reaction rates {trv) 0] for 
the first and second step, respectively. 

To calculate the second step reaction rate ^n^He), the ^IIe(n,7)^IIe cross 
section is required. One may estimate that the electric dipole transition El 
from the incident s-wave is dominating the reaction. This is confirmed by both 
three-body and direct-capture calculations. The cross section is given by 

(2) 

with 

^ Y 1 1 ^''2Ub(£^l,fl)''lMsc(£^2,?-2)f-2/fclj=3/2(?'l,?'2) , (3) 

2 i/ 

where E^ = Ei + E2 + Q12 is the energy of the photon, with Q12 being the Q- 
value of the reaction ^IIe(2n,7)^IIe. In Eq. the coordinates ri and r2 are the 
distances of the valence neutrons from the a-core and ^He, respectively. The 
quantities V2 and k2 are the relative velocity and wave number in the entrance 
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Table 3.1. Relevant input parameters for our calculations: Q12 is the Q-value for 
both steps, , Er and /r, are the spin-parity assignments, resonance energies and 
widths of the indermediate nuclei ^He, ^He and ^"Li 

reaction Q12 [MeV] j" £^r [MeV] Tr [MeV]~ 

4He(2n,7)6He O.gS'^) |"^) 0.89'^) O.TGiO.ir) 

6He(2n,7)8He 2.14'^) f"*^) 0.44^) 0.16 ±0.03*') 

^Li(2iv7)"Li 0.31^) 1+'^) 0.54'^) 0.22 ± 0.01'^) 

Rcf. 0; Rcf. ||l|; Rcf. d) Rcf. Q 



channel (^He+n) of the second step. The quantity Uh{Ei,ri)/ri represents 
the radial part of the quasi-bound wave function of the '^He-resonance, and 
Usc{E2,r2)/r2 is the radial part of the (^Hc+n)-scattcring wave function. The 
function /f=i j=3/2(?'i, ^2) is the radial component in the expansion 

^ {^Re) = ftj{ri,r2) [xtj (fi, f^i^) Xij {^2, 'T2z)]j^o (4) 



of the a+n+n ground-state wave function. Here xij sre the angular-momen- 
tum functions with orbital angular momentum i and total angular momentum 
j = £ zL 1/2, and the brackets [. . .] denote vector coupling. The expansion 
includes only equal ^-values for both neutrons and only the {£ — l,j — 3/2)- 
component contributes to the cross section. 

To calculate the first step reaction rate (n'^He) we need the (^He+n)- 
potential. It was obtained by fitting a folding-type potential |l^ to the 
experimental energy and width of the ^He ground state. For the calculation of 
the energy dependent width T (^He, Ei) of ^He in Eq. (0) we used an equivalent 
square well potential with the same r.m.s. radius of the density distribution as 
for the folding potential. 

In the second step the ground-state wave function Ub{Ei,ri) of ""He was 
determined by solving the single-particle Schrodinger equation, using the 
(^He-l-n)-potential and a boundary condition for a decaying system. It was 
then cut at a reasonable cut-off radius. For the determination of the scattering 
wave function Usc{E2,r2) we used again folding-type potentials. An accurate 
three-body wave function of ^He (see |l^) was used in our calculations. 

The relevant input parameters for the calculations of the reaction rate 
are shown in Table 3.1. The calculated cross section for ^He(n,7)®He can be 
parametrized by tr = tri /{E2 [MeV])^/^ [/ib]. The ^}ie{n,-/fRe cross section 
exhibits the well-known l/ii-behavior for incident s- waves. 

The obtained reaction rate can be parametrized in the following way: 



Nl (2n'*He) ^ aT^ exp ( ^ ) IQ-^ cm^ s"^ mol"^ , (5) 



where the parameters a, b, c are shown in Table 3.2 for the two temperature 
regions 0.1 < Tg < 2 and 2 < Tg < 15 (Tg: in units of lO^K). The inverse 
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Table 3.2. Parameters of the cross section for ^He(n,7)®He and the reaction rates for 
''He(2n,7)''He and the reverse photodisintegration of ®He calculated in our three-body 
model and in direct capture 





Three-body calculation 


Direct-capture calculation 


CTi [McV^/^ 


0.152 


0.270 




0.1 < Tg < 2 2 < Tg < 15 


0.1 < Tg < 2 


2 < Tg < 15 


a 


0.00265 


0.293 


0.00471 


0.520 


b 


2.55 


-0.351 


2.55 


-0.351 


c 


0.181 


-5.24 


0.181 


-5.24 


d 


3.14 


0.0286 


3.14 


0.0286 


e 


-1.05 


1.85 


-1.05 


1.85 


/ 


-21.9 


-16.5 


-21.9 


-16.5 



reaction rate A-y per nucleus per second can be calculated by using the RevRatio 

m 

= RevRatio x nI (2n^He) s~^ , (6) 

where we assumed that the photodisintegration proceeds through the ^He res- 
onance. The RevRatio is parametrized in the following way 

RevRatio = dT^ exp (j^j 10^^ cm"^ mol^ . (7) 



The parameters d,e, f are also shown in Table 3.2, 

The reaction rate for ^He(2n,7)^He calculated with the help of our three- 
body model (second column of Table |3.2[ ) is at least more than three orders 
of magnitude larger than the previously adopted value [|l). This is mainly 
due to the non-resonant El-transition of the second step ^He(n,7)^He, which 
dominates the cross section, and which was not taken into account previously. 



4 Direct capture (DC) calculations 

We also calculated the reaction rates for *He(2n,7)^He and the reverse photo- 
disintegration of ^Hc by using the DC model for the second step of the reaction. 
In this case Eq. ^ reduces to 

I=^^^¥^ f druUE2,ryue=,,,^y2{r) , (8) 

K2 Jq 

where r = r2 and 6*^=1^=3/2 is the spectroscopic factor and M£=i.j=3/2(?') is the 
wave function for ^He = ^He-|-n. 

The spectroscopic factor for ^He(n,7)^He was calculated by shell-model 
calculations using the (6-16) 2BME interaction of Cohen and Kurath and 
is given by 5^=1^=3/2 = 1.3139. All the other parameters (see Table |3T|) for 
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Table 4.1. Parameters for the cross section of the second step and reaction rates of 
^He(2n,7)*He and ''Li(2n,7)^^Li, and the reverse photodisintegration of *He and ^^Li 
calculated in direct capture 





SHe(2n,7)«He 


9Li(2n 


,7)"Li 


0" 1 Mcv 


1.828 


0.2435 




0.1 < Tg < 2 2 < Tg < 15 


0.1 < Tg < 2 


2 < Tg < 15 


a 


8.015 36.91 


0.651 


0.956 


b 


-0.209 -1.187 


0.042 


-0.412 


c 


-2.514 -4.262 


-3.394 


-3.432 


d 


0.013 0.003 


0.004 


0.003 


e 


1.709 2.687 


1.458 


1.912 


f 


-27.427 -25.679 


-6.470 


-6.432 



the DC calculation of '*He(2n,7)SHe were assumed to be the same as in the 
three-body calculation. For th e D C calculation we used the code TEDCA [p3| . 

As can be seen from Table |3.2| the cross section for the second step and the 
total reaction rate of ''He(2n,7)^He is enhanced by 78% enhanced compared 
to the three-body calculations. Therefore, we assumed that the DC model 
is also adequate within a factor of two to calculate the reaction rates for two 
other two-step processes involving halo nuclei like ^He(2n,7)^He and ^Li(2n,7), 
and the photodisintegration of *He and ^"'^Li. These reactions could also be of 
importance in the a-process. 

We assume as for ^He(2n,7)^He that again two p-neutrons are transferred 
in the two-step process and that in the second step the El-transition is domi- 
nating. For ^Li(2n,7)^^Li we did not consider yet possible predicted low-energy 



intruder s-states in ^"Li |24, |22l and/or ^^Li g6[. Such states could lead 
to s-wave resonances in ^^Li near the (^'^Li-|-n)-threshold modifying our calcu- 
lated non-resonant cross section for ^''Li(n,7)^^Li. Furthermore, by considering 
such s-wave states the two-step process could also proceed over 1~- or 2~- 
states in ^ "Li and not only over the l+-state at 0.54 MeV above threshold (see 
Table 3T) that we considered in our DC-calculation. 



The used input parameters are again shown in Table 3.1. The spectroscopic 
factors are given by 81^1^^3/2 — 3.2047 for ''He(n,7)^He, and S'f=i.j=3/2 = 
0.647 and S'£=i.j=i/2 = 0.024 for i°Li(n,7)"Li. For the potentials in the exit 
channels we used the folding method adjusted to the separation energy and the 
r.m.s. radii of the neutrons for the residual nuclei. 

The results for the parameters of the cross sections for the second step 



and the reaction rates are also shown in Table 4.1. The reaction rate for 



'*He(2n,7)8He is at Tg « 1 about a factor 40 larger than for 4He(2n,7)'^He. Very 
similar results for the reaction rate of '''He(2n,7)^He have also been obtained 
in another DC calculation by Gorres et al. |2^. However, our reaction rate for 
^He(2n,7)®He is more than a factor 10 larger than in ref. mainly because 
we used a larger spectroscopic factor. The reaction rate for^Li(2n,7)^^Li (Ta- 
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ble 4.1) is at Tg « 1 almost a factor two enhanced compared to He(2n,7)^He 



(Table |3.2D, even though the total Q-value for the Li-reaction is a factor three 



smaller (Table 3.1). This is mainly due to the more pronounced halo structure 



of the Li-nuclei as compared to the He-nuclei. 
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